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Single-photon superradiance from a quantum dot 
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We report on the observation of single-photon superradiance from an exciton in a semiconductor 
quantum dot. The confinement by the quantum dot is strong enough for it to mimic a two-level atom, 
yet sufficiently weak to ensure superradiance. The electrostatic interaction between the electron 
and the hole comprising the exciton gives rise to an anharmonic spectrum, which we exploit to 
prepare the superradiant quantum state deterministically with a laser pulse. We observe a five-fold 
enhancement of the oscillator strength compared to conventional quantum dots. The enhancement 
is limited by the base temperature of our cryostat and may lead to oscillator strengths above 1000 
from a single quantum emitter at optical frequencies. 


Enhancing and tailoring light-matter interaction is at 
the heart of modern quantum physics, partly because it 
enables studying hitherto unexplored realms of physics 
and partly to meet the steep requirements for quantum- 
information science. Photonic nanostructures efficiently 
tailor the density of optical states and have proven very 
useful to this end. Eor example, cavities can reach strong 
coupling to emitters [1, 2] or mechanical objects [3], and 
photonic waveguides enable efficient photonic switches [4] 
and single-photon sources [5]. Another approach to en¬ 
hancing light-matter interaction concerns tailoring the 
capability of the emitter to be polarized, i.e., the oscil¬ 
lator strength. This can be achieved with collective ef¬ 
fects such as superradiance [6], which has been studied 
in ensembles of atoms [7], ions [8], Bose-Einstein con¬ 
densates [9], and superconducting circuits [10]. Collec¬ 
tive enhancement can occur at the single-photon level if 
a single quantum of energy is distributed coherently in 
an ensemble [6]. This single-photon superradiance (SPS) 
has been studied so far in ensembles of non-interacting 
emitters such as nuclei [11], and is central to schemes 
for robust quantum communication [12] and quantum 
memories [13]. A drawback of non-interacting systems 
is their harmonic energy structure, which prohibits de¬ 
terministic preparation of a particular collective state. 
Here we show that the fundamental optical excitation of 
a weakly confining quantum dot is a generalization of 
SPS. We prepare the collective quantum state determin¬ 
istically with a laser pulse and demonstrate its superradi¬ 
ant character. Our findings underline the extraordinary 
potential of weakly confining quantum dots for achieving 
unprecedented light-matter coupling strengths at optical 
frequencies, which would improve the radiative efficiency, 
quantum efficiency, quantum nonlinearities, and coher¬ 
ence of single-photon sources in nanophotonic quantum 
devices [14]. 

We study quantum dots formed by intentional mono- 
layer fluctuations of a quantum well, which were pio¬ 
neered by Gammon et al. [15], cf. Eig. 1(a). The sub¬ 


wavelength size of the quantum dot is key to achieving a 
large collective enhancement; in larger ensembles, such as 
atomic clouds, the enhancement is reduced by destructive 
interference [6, 16]. The fundamental optical excitation 
of the quantum dot is an electron-hole pair bound by 
electrostatic attraction and quantum confinement, i.e., 
an exciton. We demonstrate that the exciton recom¬ 
bines radiatively with a quantum efficiency of (99 ± 2) %, 
which is the highest reported on quantum dots so far [17- 
19]. The resulting single photons inherit the superradiant 
character in the form of an enhanced emission rate com¬ 
pared to conventional strongly confining quantum dots. 
We employ a recently developed method exploiting the 
fine structure of the exciton [19] and measure an oscillator 
strength of up to 96 ±2. The corresponding superradiant 
enhancement of about 5 is limited by the base tempera¬ 
ture of our cryostat (7K) and could potentially be orders 
of magnitude larger at temperatures below 1 K. 

The hallmark of SPS is the symmetric collective quan¬ 
tum state [6], 

E I5i52-ej...5iv), (1) 

where N is the number of emitters, the j-th emitter is 
in the excited state |e) and all others are in the ground 
state 1^). The remarkable property of \A/) is that it in¬ 
teracts with light N times stronger than a single emitter. 
This state describes a non-interacting ensemble, where 
the excitation is localized in a single emitter at a time as 
depicted in Eig. 1(b). In a system of interacting particles, 
such as a semiconductor quantum dot, the wavefunctions 
of the underlying atoms overlap, leading to delocalized 
excitations. This destroys the collective enhancement of 
light-matter interaction and causes conventional quan¬ 
tum dots to exhibit small oscillator strengths of about 
10, despite that they embody tens of thousands of atoms. 

The spatial extent of delocalized excitations is a fun¬ 
damental property of semiconductors and is determined 
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FIG. 1: Superradiant excitons in quantum dots, (a) A quan¬ 
tum dot defined by intentional monolayer fluctuations weakly 
confine electrons (e) and holes (h), which are mutually bound 
by electrostatic attraction. Notably, the spectrum is anhar- 
monic due to interactions, i.e., the energy /kJxx of a biexciton 
is less than the energy of a single exciton. The swirling ar¬ 
row indicates superradiantly enhanced light-matter coupling, 
(b) SPS is dehned in an ensemble of non-interacting emit¬ 
ters as a symmetric superposition of different excitations, (c) 
The excitonic enhancement of light-matter interaction may 
be regarded as a generalization of SPS: the exciton is a sym¬ 
metric superposition of excitations, (d) Measured photolu¬ 
minescence spectrum at 10% of the exciton saturation power 
Psat = 20 nW. Only the exciton is observed, (e) At the sat¬ 
uration power, the biexciton becomes discernible, (f) The 
excitons and biexcitons are distinguished by their power-law 
dependence on excitation power, P\ the fits yield and 

respectively. 


by the size of an exciton. Enhancement of light-matter 
interaction can therefore be achieved only in quantum 
dots that are larger than the exciton radius. This regime 
is known as weak confinement and the enhancement of 
the light-matter coupling in weakly confining quantum 
dots was first predicted by Hanamura [20] and dates back 
to early theoretical studies of impurities in semiconduc¬ 
tors [21]. Here we show that this effect is equivalent to 
SPS, and the exciton state can be written as, see Supple¬ 


mentary information [22], 
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where r (R) is the relative (center-of-mass) electron-hole 
coordinate, and the index j runs over the unit cells con¬ 
stituting the quantum dot. The function 0x describes an 
exciton with the size and oscillator strength of a conven¬ 
tional quantum dot, while c is responsible for the collec¬ 
tive enhancement as illustrated in Fig. 1(c). The light- 
matter coupling is proportional to the number of atoms 
comprising the weakly confining quantum dot. This anal¬ 
ysis shows that this effect is a generalization of SPS, cf. 
Figs. l(b,c), and the two effects are equivalent if c is con¬ 
stant throughout the quantum dot. The constant phase 
of c found for ground-state excitons with 5-like symmetry 
ensures constructive interference among the excitations 
defined by 0x- 

An exciton governed by Eq. (2) has been long sought 
in solid-state quantum optics [23] because it can lead 
to large oscillator strengths. Realizing weakly confining 
quantum dots has been a challenge so far because it re¬ 
quires precise control over growth parameters to obtain 
a homogeneous potential profile over extended length 
scales. Previous studies on large quantum dots [18, 19] 
revealed small oscillator strengths, which is believed to 
be caused by inhomogeneous potential profiles within the 
quantum dots. Previous works [24, 25] found fast re¬ 
combination rates in gallium-arsenide quantum dots but 
without rigorous information about the impact of non- 
radiative processes, see Supplementary information [22] 
for further discussions. The measured spectra shown in 
Fig. 1(d) and (e) were obtained by exciting in the quasi¬ 
continuum energy band of the quantum well as discussed 
below. An exciton and a biexciton are identified as shown 
in Fig. 1(f). These quasi-particles radiate at different fre¬ 
quencies, cf. Fig. 1(e), which reflects the spectral anhar- 
monicity of the quantum dot. 

To identify proper excitation conditions of the quan¬ 
tum dot, we probe the spectrum of states with 
photoluminescence-excitation spectroscopy as displayed 
in Fig. 2(a). The spectrum shows a quasi-continuum 
band of quantum-dot states hybridized with quantum- 
well states as well as the exciton manifold in which we 
identify the Is, 2s, and 3s states of two-dimensional ex¬ 
citonic hydrogen [26]. Key features of the spectrum are 
summarized in Fig. 2(b). We use two excitation condi¬ 
tions to prepare the Is exciton: i) Pumping in the quasi¬ 
continuum band of states (C-type excitation) allows ex¬ 
tracting the impact of non-radiative processes govern¬ 
ing the Is exciton decay. Since the quantum dot traps 
carriers with random spin, equal populations of spin- 
bright and spin-dark Is excitons are prepared. While 
only bright excitons emit light, the dark excitons in¬ 
fluence the decay dynamics and play a key role in re¬ 
vealing non-radiative effects (see Supplementary infor- 
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FIG. 2: Deterministic preparation of superradiant excitons. 
(a) Photoluminescence-excitation spectrum obtained by inte¬ 
grating the emission of the Is transition while scanning the 
excitation wavelength. It features a quasi-continuum and re¬ 
solves the lowest-energy states of the exciton manifold, la¬ 
beled Is, 2s, and 3s. (b) Two excitation schemes are used 
in our study. With C-type excitation, an equal bright- and 
dark-exciton population is prepared, which is important for 
extracting the impact of non-radiative processes. Determinis¬ 
tic preparation of the bright exciton is achieved by pumping 
into the 2s state. 

mation [ 22 ]). ii) Deterministic preparation of spin-bright 
superradiant Is excitons is achieved by pumping into the 
2s exciton state, cf. Fig. 2, with a pulsed laser. This 
is feasible since the decay cascade from 2 s to Is is spin- 
conserving [27] and spin-dark states are not populated. 
Deterministic excitation occurs when applying sufficient 
optical power (300 nW) to saturate the emission from the 
Is state. 

The figure of merit for collective enhancement of 
light-matter interaction is the oscillator strength, /, 
which gauges the strength of the interaction with light. 
The oscillator strength is determined by the radiative 
spontaneous-emission rate of an emitter placed in a ho¬ 
mogenous photonic environment. In an experiment, how¬ 
ever, the oscillator strength is masked by non-radiative 
effects and the non-homogeneity of the photonic envi¬ 
ronment, whose contributions are fully addressed in our 
study, see Supplementary information [ 22 ] for further de¬ 
tails. Central to our analysis is a recently developed 
method exploiting the fine structure of excitons to rigor¬ 
ously separate radiative from non-radiative effects [28]. 
Figure 3a shows raw data of the time-resolved decay of 
the deterministically prepared Is exciton. We measure an 
excellent near-unity radiative efficiency of 77 = (99±2)%, 
which is the highest ever measured on quantum dots so 
far. The extracted oscillator strength of / = 72.0 ± 0.8 
is enhanced far beyond the upper limit of / = 17.4 for 
conventional quantum dots at this wavelength. By com¬ 
bining structural information about the sample with the 
measured oscillator strength we can faithfully reconstruct 
the exciton wave function and find a diameter of 24 nm, 
which is smaller than the wavelength of light yet suffi¬ 


ciently large to embody ^ 90, 000 atoms in a collective 
quantum state sharing a single quantum of energy. 

Figure 3(b) shows the second-order correlation func¬ 
tion obtained in a Hanbury-Brown-Twiss (HBT) exper¬ 
iment from which we find a normalized zero-time corre¬ 
lation function of ^^^^( 0 ) = 0.13. Antibunched emission 
directly pinpoints to the presence of a single quantum 
of matter inside the quantum dot, which is an essential 
property of SPS, cf. Eq. ( 1 ). In conjunction with the mea¬ 
sured enhanced oscillator strength for a spatially confined 
exciton, this is the unequivocal demonstration of SPS in 
a quantum dot. Solid-state quantum light sources often 
suffer from charge traps that switch the emitter into an 
optically dark state, also known as blinking processes, 
and reduce the preparation efficiency of bright states. 
This can be quantified from HBT correlations acquired 
over long time scales as shown in Fig. 3(c). No bunching 
effects are observed, which shows that this single-photon 
source is free from blinking on a time scale of at least 
10 ps. 

Single-photon superradiance is a robust phenomenon 
in quantum dots due to the anharmonic spectrum. The 
only experimentally relevant parameter that may be 
detrimental to SPS is temperature. For a thermal de 
Broglie wavelength larger than the quantum-dot size, 
the exciton Pauli blockage is broken and leads to multi¬ 
photon emission, thereby destroying SPS. This non¬ 
trivial effect is beyond the scope of the current work and 
will be presented elsewhere. 

We have measured the oscillator strength of 9 quan¬ 
tum dots and found them all to be superradiant with 
an average oscillator strength of / = 76 ± 11. Fur¬ 
ther experimental data is included in the Supplemen¬ 
tary information [ 22 ]. Remarkably, we have measured 
a homogeneous-medium radiative decay rate of up to 
Prad = ( 11.1 ± 0 . 2 )ns“^, which is the fastest value ever 
reported for any single-photon source and corresponds to 
an oscillator strength of / = 96 ± 2 . Such a quantum dot 
can deliver a radiative flux of single photons equivalent 
to more than five conventional quantum dots. 

The superradiant enhancement of the light-matter cou¬ 
pling in quantum dots is proportional to the number of 
atoms in the collective state, and can potentially be much 
larger than reported here. The enhancement factor may 
realistically reach 100 x for quantum-dot diameters of 
^ 100 nm [29] corresponding to an oscillator strength of 
/ ^ 1500. Such highly superradiant quantum dots may 
exist in our sample but the temperature at which the 
experiment is carried out (T = 7K) does not allow re¬ 
solving such large oscillator strengths. This is because in 
large quantum dots the confinement energy may become 
smaller than the thermal energy, which results in pop¬ 
ulating excited states with reduced oscillator strength. 
The maximum oscillator strength /max,th that can be re¬ 
solved at a temperature T is calculated for a quantum 
dot in which the energy difference between two eigen- 
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FIG. 3: Experimental demonstration of single-photon superradiance from a quantum dot. (a) Time-resolved decay (black 
points) of the bright Is exciton obtained under 2s excitation. We obtain an excellent fit (yellow line) to the theoretical model 
when convoluting with the instrument-response function of the detector. After separating non-radiative from radiative effects 
we extract a radiative decay rate of (8.3 ± 0.1) ns“^ (red line), which is deeply in the superradiant regime (green area), (b) 
HBT measurement of the emitted photons showing ^'^^^(0) = 0.13, which proves single-photon emission, (c) Long-time-scale 
HBT measurement where each coincidence peak has been numerically integrated. No blinking is observed. 


states equals d/cpT (a detailed analysis is provided in the 
Supplementary information [22]) 


_ AhEp 1 1 


(3) 


where Ep is the Kane energy, M the exciton mass, ao 
the exciton radius, kpT the thermal energy, and 1:^ 
the in-plane aspect ratio of the quantum dot (^ > 1). 
At T = 7K we find that oscillator strengths larger 
than /max,th = 170 cannot be resolved for in-plane 
symmetric quantum dots, and /max,th decreases even 
further for more realistic asymmetric shapes. Oscil¬ 
lator strengths of ^ 1500 require temperatures below 
^ 0.8 K. The light-matter coupling may be further en¬ 
hanced > lOx beyond the homogeneous-medium value 
by the Purcell effect [14]. This could allow studying fas¬ 
cinating non-energy-conserving effects such as the ultra¬ 
strong regime of light-matter coupling [30]. The repeti¬ 
tion rates of single-photon sources would approach the 
terahertz regime yielding radiating powers of hundreds 
of nanowatts from a single quantum emitter. The single¬ 
photon emission would potentially be highly coherent, 
partly due to an intrinsically weaker coupling to nuclear 
spin noise [31] and phonon dephasing [32] for large exci- 
tons, partly because the dephasing mechanisms present 
in solid-state environments would be largely negligible 
compared to a radiative decay at subpicosecond time 
scales. Even larger decay rates could become possible 
in materials with small Bohr radii [33, 34], see Eq. (33). 
In particular, fast decays have recently been reported in 
CdSe nanoplatelets [35] . Another intriguing aspect of the 


SPS regime is that the collective Lamb shift is predicted 
to be finite [36] without the renormalization schemes re¬ 
quired in the quantum electrodynamics of conventional 
emitters. 
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I. THEORY OF SINGLE-PHOTON SUPERRADIANCE FROM QUANTUM DOTS 


Quantum dots are well suited for enhancing the light-matter interaction strength due to their multi-atomic nature. 
However, quantum dots that are smaller than the exciton Bohr radius exhibit a limited oscillator strength because 
the electron-hole motion is uncorrelated and dominated by quantum-confinement effects. This is known as the strong- 
confinement regime. In the opposite limit, weak confinement, the relative electron-hole motion is strongly correlated 
while their center-of-mass motion is bound to the weakly confining potential of the quantum dot [ 20 ]. 


The strong-confinement regime 

In the strong-confinement regime, the exciton wave function can be written as a product of the individual wave 
functions of the electron and the hole, i.e., the electron-hole motion is decoupled as illustrated in Fig. 4(a). In this 
limit the oscillator strength is given by [19] 


f=j^m{r)\F,{r))\\ (4) 

where Ep is the Kane energy and F’e(r) and Fh{r) are the slowly varying envelopes satisfying the single-particle 
effective-mass Schrodinger equation for the electron and the hole, respectively. The oscillator strength of small 
quantum dots has therefore an upper limit of /max = Epjhuj amounting to 17.4 for GaAs quantum dots at a wavelength 
of 750 nm. 


The weak-confinement regime 

Excitons in large quantum dots are bound by the mutual electrostatic attraction between electrons and holes. In 
this regime, which applies when the mean electron-hole distance, i.e., the exciton Bohr radius, is smaller than the 
quantum-dot size, the exciton wave function is non-separable implying that the electron and the hole are spatially 
entangled. The oscillator strength is proportional to the volume of the quantum dot. This effect is sometimes referred 
to as the giant oscillator-strength effect for quantum dots (GOSQD) and was first considered by Hanamura [20]. 
Excitonic enhancements of the light-matter coupling has been studied theoretically in a range of solid-state systems 
including impurities in bulk semiconductors [ 21 ] and quantum wells [37]. The unique feature of the GOSQD is that 
it occurs for a discrete quantum state. 

Our experiments concern investigations of single GaAs interface-fluctuation quantum dots embedded in 
Alo. 33 Gao. 67 As as presented in Eig. la of the main text. Bound excitonic states are obtained by intentionally cre¬ 
ated monolayer fluctuations in a quantum well, leading to a weak in-plane quantum confinement. The quantum-well 
thickness is smaller than the bulk exciton Bohr radius of ao = 11.2 nm [38] and leads to strong confinement in the 
growth direction. Exciton enhancement is achieved only within the plane, where the quantum-dot wave function is 
extended beyond the exciton Bohr radius. We model this as a cylindrically symmetric quantum dot whose slowly 
varying envelope is separable into in-plane and out-of-plane components and the out-of-plane component is further 
separable in independent components for the electron and the hole. Hence, the electron-hole motion is correlated in 
the plane and uncorrelated perpendicularly and the total excitonic wave function can be written in the effective-mass 
approximation as 


4^(R, r, re, r/,) = T)2pe{ze)iJh{zh)ue{rh)ux{re), (5) 

where i*) is the in-plane slowly varying envelope, ' 0 e(^e) (' 0 e(^/i)) is the out-of-plane envelope function for the 

electron (hole), Ue (ux) the electron (heavy-hole) Bloch function at the F point in reciprocal space, R = {rrieVe + 
+ ^/i): I* = Tg — thc ccntcr-of-mass and relative in-plane coordinates of the exciton, and rrie and ruh 
are the electron and hole effective masses, respectively. The unit-cell Bloch functions contribute to the Kane energy 
and do not play an important role in our study, which is why only the slowly varying component 4^x is addressed in 
the main text. The slowly varying envelopes in the growth direction can be accurately computed because the 
quantum-dot thickness is known precisely and amounts to Lz = 4.3 nm but they play no role for the GOSQD effect, 
which is governed by the in-plane excitonic envelope 4 ^x(H 5 i*). To see this, we first make some realistic assumptions 
and consider symmetric in-plane parabolic quantum confinement, in which case the excitonic envelope separates into 


9 




FIG. 4: Superradiance with single quantum dots, (a) In small quantum dots, such as self-assembled In(Ga)As quantum dots, 
the motion of electrons and holes is governed by quantum-confinement effects and is completely uncorrelated, which limits the 
light-matter interaction strength, (b) In large interface-fluctuation quantum dots, the electron-hole motion is dominated by 
their mutual attraction in the plane of the quantum dot. The mean electron-hole separation exciton Bohr radius; dark gray) 
is smaller than the exciton wave function (green) and leads to a strong superradiant behavior of the ground-state exciton. 
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FIG. 5: Galculated superradiant oscillator strength for an interface-fluctuation quantum dot normalized to the strong- 
confinement limit of /max = 17.4. 


center-of-mass xcm(R) and a relative-motion Xr(i*) components [29, 39] 


5'x(R.,r) = xcM(R.)Xr(r) 

(6) 

:,cm(R) = 

V TT yO 

(7) 

x.(r) = .A J-e-IH/^Qw, 

V TT agw 

(8) 


where agw is the exciton Bohr radius in the quantum well and (3 the in-plane size of the quantum dot. For a perfect 
two-dimensional system, the exciton Bohr radius is twice as small agw = tto/2 5.6 nm leading to a binding energy 
four times as large. The structure investigated in this work is, however, not a perfect two-dimensional system because 
the exciton wave function has a non-zero thickness. As argued in Refs. 40, 41, the binding energy of an exciton in a 
4 nm-thick quantum well is only twice larger than in bulk. We therefore consider a value of the two-dimensional Bohr 
radius agw — ao/\/2 « 8nm. For P > the mean distance between the electron-hole pair (« 2agw) is smaller 

than the quantum dot size (~ 2/3). 


Relation between the GOSQD effect and single-photon Dicke superradiance 

The connection to the single-photon Dicke superradiance can be made by noting that, \i ^ > aqw, the center-of- 
mass motion can be written as a convolution between a function Ca(R) capturing the dynamics on the scale of agw 
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and a function Cs(R) responsible for the coherent superradiant enhancement, i.e., 

Xcm(R.) =c„(R)*Cs(R) = /'d2pcs(P)c„(R-P) 

(9) 

«^c(Rj)Ca(R-Rj), 

j 

where the last step involves switching the integral to a sum over unit cells and c equals Cg times the discretization 
area. Consequently, the slowly varying excitonic envelope reads 


4'x(R,r) = y]c(Rj)</>x(R--R.i,r), (10) 

3 

where 0x(R, i*) = Ca(R)Xr(i*) is the two-dimensional exciton wave function and is identical to from Eq. (2) in 
the main text. The physical meaning of Eq. (10) is illustrated in Eig. lb in the main text and the analogy with the 
Dicke superradiance becomes clear, if compared with Eq. (1) in the main text. Eor parabolic in-plane confinement we 
obtain the following expressions for and Cg 


Ca(R) = 2—6 2 |R|Voqw 

(11) 

7rag 

Cs(R) = , 

(12) 


where ^ for P uqw Since the phase of Cg is constant throughout the quantum dot, the ground- 

state exciton is found in a spatial superposition with constructive cooperativity and enhanced coupling to the light 
field. In the following we quantify the expected superradiant increase in the oscillator strength and, consequently, in 
the spontaneous-emission rate. 

According to Eermi’s Golden Rule, the probability of photon emission is proportional to | (0 \px \ |^- 

The relative motion is taken to be zero, r = 0, because the exciton can recombine radiatively only if the electron and 
hole are found at the same spatial position [29]. After performing the standard procedure of merging the unit-cell 
Bloch functions into the Bloch matrix element Pcv = E[7c \Px \ '^e)uc’ where the subscript UC denotes integration 
over a unit cell, we obtain the following expression for the oscillator strength (compare with Eq. (4)) 

/ = ^X.(O) \{0\xcM{R))f \{Mz)mz))f , (13) 

where the first (second) inner product denotes a two-dimensional (one-dimensional) integration over R {z). We define 
the radius of the quantum dot L = ^/2/3 as argued in Ref. 18 and, with the help of Eqs. ( 6 - 8 ), arrive at the following 
superradiant enhancement S of the oscillator strength 

^ = 7^ = i^] ( 14 ) 

/max J 

Note that this equation is valid only in the weak-confinement regime where [3 ^ < 2 qw- The electron and hole 
wave functions in the growth direction can be accurately calculated and we find that |('0/i|'0e)| ~ 0.96 for the 

interface-fluctuation quantum dots from the present study. We plot the resulting superradiant enhancement of the 
oscillator strength in Eig. 5. It scales with the quantum-dot area and is a dramatic effect; for realistic quantum dot 
diameters of 35 nm, the light-matter interaction strength exceeds the upper limit of strongly confined quantum dots 
with uncorrelated electron-hole pairs by an order of magnitude. 


II. SAMPLE AND EXPERIMENTAL SETUP 

The sample used in our experiment was grown on a GaAs (001) wafer following the procedure developed by Gammon 
et al. [15]. The GaAs interface-fluctuation quantum dots were created by random monolayer fluctuations in the GaAs 
quantum-well thickness. The GaAs quantum well is surrounded by 5-nm-thick AIq. 33 Gao. 66 As layers in order to obtain 
a high-quality interface, and followed by a 100 -nm-thick AIq. gGao .2 As. The detailed structure is presented in Eig. 6 (a). 
A zirconia solid-immersion lens shaped as half a sphere with a radius of 1 mm and refractive index of 2.18 was placed 
on top of the sample to improve the collection efficiency. 




11 



1 mm 

10nm 
100 nm 
5 nm 

(3.7 ± 0.3) nm 
5 nm 

lOOnm 


(b) 

For power 
stabilization 


Detection 



PM fiber 
Polarizers 


PM fiber 


Photodiode J 






it 



Cryostat 


Excitation 

Vaccum window 

Microscope 
-objective 

— Heater 

Piezo-stages 


(c) 


Flip mirror CCD 






FIG. 6: (a) Cutaway profile of the investigated sample (not to scale). Lattice-matched GaAs quantum dots are formed by 

random fluctuations in the GaAs quantum-well thickness. A zirconia solid-immersion lens enhances the collection efficiency 
of the setup, (b) Sketch of the optical setup around the cryostat. Optical excitation and collection are performed in a cross- 
polarized scheme to discriminate between the photoluminescence and specular laser reflection, (c) The emission is sent through 
a spectrometer before being detected by a GGD for spectral measurements and an APD for time-resolved measurements, (d) 
For correlation measurements, the quantum dot emission line is spectrally filtered before being directed onto a Hanbury-Brown- 
Twiss (HBT) setup. 


There are several types of optical measurements performed in this study: spectral and time-resolved measurements, 
and second-order correlation measurements. All of them are carried out in a closed-cycle cryogen-free cryostat as 
sketched in Fig. 6(b). The sample holder is mounted on piezoelectric nanopositioning translation stages. For all 
experiments, the sample is cooled to a temperature of 7K. After exiting the single-mode polarization-maintaining 
(PM) fiber, the excitation beam generated by a picosecond pulsed TiiSapph laser is collimated to a diameter of 2 mm. 
Then, it passes through a thin-film linear polarizer and a 90:10 (transmitted:reflected) beam splitter before being 
focused on the sample through a microscope objective with a numerical aperture of 0.85. The spatial resolution of 
the objective was measured to be 1.1 pm^ at a wavelength of 633 nm. The excitation laser is tuned to a wavelength 
of about 750 nm corresponding either to resonant excitation of continuum states in the quantum well or to 2s-shell 
excitation of the quantum dot. The photoluminescence of the investigated ground-state excitons is located around 
752 nm. The emission is collimated by the same microscope objective and filtered from the excitation laser by the 
perpendicularly-oriented thin-film linear polarizer, see Fig. 6(b). The beam is then coupled into a PM fiber and guided 
towards the detection setup. 

Spectral measurements are performed by sending the emission to a spectrometer with a groove density and spectral 
resolution of 1200mm“^ and 25pm, respectively, and subsequently detected by a charge-coupled device (CCD), see 
Fig. 6(c). After the grating, a mirror can be flipped to direct the emission to an avalanche photo-diode (APD) with 
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FIG. 7: Comparison between photoluminescence (PL) and photoluminescence excitation (PLE) spectra from QD2. (a) PL 
spectrum obtained from exciting the light-hole resonance in the quantum well at 80 nW. The red-shaded area corresponds to the 
data integrated to obtain the PLE. (b) PLE spectrum performed across the main quantum dot and quantum-well resonances 
at 400 nW. 


a time resolution of 60 ps for time-resolved measurements. For correlation measurements, the emission is first filtered 
by a grating with a groove density of 1200 mm“^ before being coupled back into a single-mode PM fiber and directed 
towards a beam splitter, see Fig. 6(d). The grating setup has a spectral resolution of 50 pm. After the beam splitter, 
two APDs detect coincident counts. 


III. EXCITATION SCHEMES 

In the present work we have studied 9 different interface-fluctuation quantum dots, which are labelled as 
QD1...QD9. Most of the results presented here stem from QDl except Fig. 7 (QD2), and Fig. 14 (QD3). All 
data in the main text stem from QDl except Fig. 2, which was obtained from QD2. We note that all nine quantum 
dots had very similar properties, cf. Table 1. 

To acquire a thorough understanding of the level structure of the quantum dot, we measure the electronic density 
of states in a photoluminescence-excitation (PLE) experiment. The measurement was carried out in continuous-wave 
mode below the saturation power of the Is exciton. The laser was scanned stepwise from around 752 nm down to 
735 nm where the quantum-dot and quantum-well resonances are present. For each acquired spectrum, the exciton 
photoluminescence (PL) was fitted by a Lorentzian function and integrated, and the resulting quantity is shown by the 
height of the quantum-dot line at the given spectral position of the laser as shown in Fig. 7 (b). We observe a couple 
of quantum-dot resonances stemming from the 2s- and 3s-shells before the onset of a spectrally continuous absorption 
of the quantum well. The wavelength used for quasi-continuum excitation is denoted by the C-line in Fig. 7 as well as 
in the main text. The quantum-well resonances consist of 2D continuous exciton states and quasi-continuum states 
























13 


[24]. 

The oscillator strength characterizes the coupling strength of a two-level system to light. The quantum dot excitation 
scheme needs, therefore, to be as clean as possible, such that the environment of the quantum dot is not polluted by 
phonons and charges that may raise the effective temperature at the quantum dot position and couple the excitonic 
levels thermally. This requirement is particularly stringent for interface-fluctuation quantum dots owing to the close 
proximity of the resonances of only a few meV induced by weak quantum confinement. We therefore use the 2s-shell 
as the main excitation scheme, which also allows us to deterministically prepare the bright superradiant state as 
explained in the main text. We note that the spectrum remains clean also when pumping through the quantum-well 
quasi-continuum, cf. Fig. l(c-d) in the main text, but the radiative decay rate is found to depend on excitation power, 
which is attributed to the presence of undesired thermal processes induced by the local phonon bath created by 
the relaxation of the (many) charge carriers from the quantum well. This excitation scheme is therefore only used 
for extracting the decay rate of nonradiative processes happening in the quantum dot as discussed below, while the 
oscillator strength is probed through the 2s excitation. 


IV. MEASURING THE OSCILLATOR STRENGTH AND PREVIOUS WORK ON GOSQD 


Experimentally, the oscillator strength can be probed either in absorption [24] by extracting the polarizability or in 
emission by extracting the homogeneous-medium radiative decay rate [19]. The former approach has the advantage 
that it is unaffected by non-radiative processes but a main drawback is that the quantized character of the excitation 
is not probed. We therefore study the radiative decay, where the quantized character is encoded in the statistics of 
the emitted light. The radiative decay is potentially masked by non-radiative and spin-flip processes [17], which are 
ubiquitous in a solid-state environment. The issue is that a single-exponential decay curve contains contributions from 
both radiative and non-radiative processes whose sum is measured by fitting the decay curve. Separating radiative 
from non-radiative processes in such an experiment is therefore not possible. In the present study we address this 
issue and carefully map radiative and non-radiative contributions using the fine-structure splitting of quantum dots, 
a technique that results in a bi-exponential decay dynamics of the quantum-dot exciton and is described thoroughly 
in Refs. [14, 19, 42]. 

The decay rate measured in an emission experiment is the total rate. 


r(r,u;,ep)=r|^,T(c.) 


p(r,cj,ep) 
Phom (^) 


■r nrad? 


(15) 


where p(r, cj, Gp) and Phom(<^) are the local density of optical states in the given photonic structure and a homogeneous 
medium, respectively, rnrad(^^) is the rate of non-radiative recombination, hw the emission energy, r the position of 
the quantum dot, and rj^°^(cd) the radiative decay rate of a quantum dot in a homogenous medium. The oscillator 
strength is given by [18] 


/(^) = 


Girmoepcl 


■phom 

rad 


M, 


(16) 


where n{(jj) is the refractive index of the material surrounding the quantum emitter, uj and Cp the frequency and 
speed of light, respectively, eo the vacuum permittivity, mp the electron mass, and q the elementary charge. We stress 
that the oscillator strength is proportional to rjl°™(u;), which can be very different from the measured decay rate 
r(r,cd,ep). We have precise knowledge of the structure of our sample (see section ) and, with the help of Ref. 43, 
calculate a small radiative inhibition of p(r, cd, ep)/phom = 0.95 compared to a homogeneous medium. The fraction 
of decay events resulting in photon emission is the quantum efficiency defined as 


7]{UJ) 


phom 
^ rad 


(w) 


^rad”(^) rnrad(^) 


(17) 


Previous searches for the GOSQD effect was inspired by the prediction by Andreani et al. [23] that quantum dots 
in the GOSQD regime may enable reaching the strong-coupling regime of cavity quantum electrodynamics. In some 
works [2, 44], the oscillator strength was estimated from the vacuum Rabi splitting in the strong coupling regime 
of cavity quantum electrodynamics. Such estimates are inaccurate because multiple quantum dots may couple to 
the cavity even when they are off resonance due to a (non-Dicke, non-single-photon) collective coupling of multiple 
quantum dots to the cavity mediated by phonon coupling [45, 46]. In other works, the oscillator strength was estimated 
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from absorption experiments [24] but in such experiments the influence of other emitters cannot be ruled out. The 
oscillator strength has also been estimated from time-resolved measurements [25, 47] but, as pointed out above and 
also noted in Ref. 25, the non-radiative and radiative processes must be measured independently. It is also crucial 
to extract properly the radiative decay rate for a homogeneous medium because the local density of optical states is 
modified significantly in photonic nanostructures even by the presence of nearby planar surfaces. The importance of 
properly accounting for these effects was highlighted in recent results on large InGaAs quantum dots: in Ref. 47, the 
total decay rate was used to estimate an oscillator strength of ^50 but later measurements showed that non-radiative 
processes were very significant and that the oscillator strength was ^5 times smaller [18], i.e., below the GOSQD 
regime. 

In order to extract accurately the oscillator strength of the quantum dots, we use time-resolved spectroscopy on 
single emitters along with an appropriate model for the bright exciton decay [18, 19]. Excitons confined in quantum 
dots arise from the binding of an electron from the conduction band with a heavy-hole from the valence band because 
compressive strain and confinement shifts the light-hole band to higher energies [48]. The electron-heavy hole complex 
has total angular momentum j = 2 with four possible combinations: rrij = {±2,±1}. Excitons with low projected 
momentum mj = ±1 are called bright \b) since they couple to light, while mj = ±2 states are dark \d) and do not 
couple to light. There is a splitting in energy Abd between the two excited states of several hundred peV due to 
electron-hole exchange interaction [49]. The level scheme of the exciton is pictured in Eig. 8, where \g) denotes the 
ground exciton state. The bright exciton can decay both radiatively and non-radiatively while the dark exciton can 
only decay non-radiatively. The exciton can flip its spin with the rate Tgf. Spin-flip processes are phonon-mediated, 
and since ksT ^ Abd for all measurements, the spin-flip rate is assumed to be the same either way between bright 
and dark excitons. This way, a dark state can contribute to the radiative decay of the bright exciton only if it first 
undergoes a spin flip to the bright state, thereby providing the slow rate. Einally, the non-radiative rates are taken 
to be the same for both excitons due to the small energy splitting between the two states [42] . 



FIG. 8: Three-level exciton scheme with the ground state \g) corresponding to no exciton, and two excited states: bright state 
\b) and dark state \d) split in energy and coupled through spin-flip interaction with the rate Fgf. Only \b) can decay radiatively 
while both states can also decay non-radiatively. 

The population probabilities of bright pb and dark pd excitons are governed by the following system of equations[17] 

= —(brad + Enrad + rsf)pb + TgfPd, (18) 

= “(Enrad + rsf)pd + TgfPb, (19) 

and is solved for the bright state pb 

Pb(t) = ('Q)g- (rnrad+FsOt ^ |^20) 

J- rad 

where Pb(0) and Pd(0) are the initial populations of dark and bright excitons, respectively. Erom time-resolved 
spectroscopy where the decay of the bright exciton is probed, we can retrieve the three decay rates E^ad, Enrad and 
Egf through 


Erad = Ef - Eg 

Pb(0) 


- nrad 


= r«- 


Tsf = 


Pd( 0 ) 

As /9b (0) 


(rf-b) 


A( Pd(0) 


(Ff - Ts 


( 21 ) 

( 22 ) 


(23) 
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FIG. 9: Decay curve of QDl taken under quasi-continuum excitation (blue) along with the bi-exponential fit (red) and the IRF 
(green). Below are the normalized residuals yielding = 1.13. The extracted parameters for the fit are indicated. 


If the initial preparation probabilities of the bright and dark excitons are known, the radiative rate can be unam¬ 
biguously extracted from the experimental data [19]. These probabilities are equal when pumping in the quantum 
well pb(0)/Pd(0) — 1 because carriers with random spins are captured by the quantum dot. If for some reason (e.g., 
spin-conserving cascade to the ground-state exciton) pb(0)/Pd(0) > 1, we would actually be overestimating Fnrad be¬ 
cause of Eq. (22). This means that we are estimating a lower bound to the oscillator strength and quantum efficiency. 
Quantum-well excitation yields the fast rate F^^ = Fj?^^ + Enrad + Egf, where the three rates correspond to radiative 
decay, nonradiative decay and bright-to-dark-state spin flip, respectively [19]. We extract Fnrad = (130 ± 3) ps“^ and 
Esf = (360 ± 30) ps“^, quantities that are not expected to depend on excitation conditions, and this is confirmed 
experimentally by measuring an excitation-independent slow rate. We note that, due to the presence of thermal 
processes, is not related to the oscillator strength but rather to an effective transition strength [19] because more 
than two levels are involved in the exciton dynamics. 

The oscillator strength is measured using 2s-shell excitation, where the fast rate is F^^ = Frad + Enrad + Egf and we 
obtain Erad = (8-3 ± 0.1) ns“^, which is more than four times as fast as the uncorrelated limit of 1.93 ns“^, see Fig. 3 
in the main text. By taking into account the optical density of states in our structure, this yields a homogeneous- 
medium decay rate of 8.7 ns“^. The experimental proof of SPS is completed in a Hanbury-Brown-Twiss experiment, 
where the signal is sent to a beam splitter and subsequently detected by two APDs located in the transmission ports 
of the beam splitter. The resulting coincidence counts on the APDs are shown in Fig. 3 in the main text. We 
note that 2s-excitation changes the initial preparation condition Pb(0) ^ 1 and Pd(0) ^ 0 owing to the 2s-ls spin- 
conserving cascade, which allows deterministic preparation of SPS. This is confirmed in the experimental data, where 
the quantum-dot decay curves under 2s excitation are reliably fitted with a single-exponential model except in a few 
quantum dots where the spin-flip rate is fast enough to intermediately create a very small dark-exciton population. 

Finally, we discuss the procedure for fitting the decay curves. In time-resolved measurements, the instrument 
response function (IRE) has to be taken into account when fitting the data. The measured signal /meas(Q is the 
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convolution of the input signal f{t) with the IRF r{t) of the APD 

t 

/meas(i) = j r{T)f{t - t)At. 

— OO 


(24) 


The laser pulse has a length of 3 ps, and is much shorter than the APD resolution. Therefore, the IRF is measured by 
collecting light from the laser reflection from the sample surface. Finally, the data are fitted by the aforementioned 
model (convoluted with the IRF), described by a sum of Ne (either 1 or 2) exponents added to a background value 




^fit(^) — BG + E hie 


-Tit 


1=1 


(25) 


where BG is the background level determined by the background counts measured on the APD, and the after-pulsing 
probability function of the count rate and wavelength. In order to estimate accurately the background level for each 
measured decay curve, we calibrated the after-pulsing probability close to the quantum dots wavelength for different 
count rates. The decay curves are fitted by a least-squares approach from which a weighted residual is estimated for 
each data point /c. 




3as (^/c) -hit(^/c) 

\/-^meas (^/c) 


(26) 


where /meas(^) is the measured data and /fit(t) is the fitted value. Finally, the sum of the squared residuals = 
C^k=i ~ P) i^ minimized to render the best fit to the data, where N is the number of time bins and p 

is the number of parameters in the model. Figure 9 shows the decay dynamics excited through the quantum-well 
quasi-continuum and the results of the fitting. 


V. UNCERTAINTY OF THE EXTRACTED PARAMETERS 


We have carefully addressed the experimental uncertainties for the quantum dot presented in the main text, see 
also Sec. . We have performed the experiment four times in different days and have done a statistical analysis 
on the extracted quantities, which accounts for experimental errors such as laser-power fluctuations, temperature 
fluctuations, etc., as well as for the sanity of the fitting routine by randomizing the error signal. The errors of the 
extracted parameters (decay rates and amplitudes) have been combined in quadrature and propagated using Eqs. (21- 
23) in textbook fashion [51]. We obtain an uncertainty of about 2% for the oscillator strength and superradiant 
enhancement as well as for the quantum efficiency. 

For the other eight quantum dots a single time-resolved measurement was performed and the uncertainties have 
been calculated by performing a thorough error analysis on the experimental data. The resulting errors are listed in 
Sec. . In the following we present the method used to extract these uncertainties. 

The bi-exponential function from Eq. (25) containing 4 parameters p = (Fj, F^, 6/, 6^)^ (note that the background 
signal BG is a fixed parameter determined by the dark counts and afterpulsing of the APD) is fitted to the raw data 
with a nonlinear regression employing the Levenberg-Marquardt iterative linearization method [52]. The regression 
finds the optimal set of values po that minimizes the sum of the squared residuals X^in- estimate the accuracy of 
Po, the variance-covariance matrix coy{p} is evaluated at po using the technique presented in Ref. 53. The square root 
of the diagonal components of this matrix denote the fitting uncertainty of p around po, i.e.. 


O'(^lp=po = \/diag[cov(^] 


(27) 


where a denotes the standard deviation. The entries in po may have different physical units and it is therefore 
more convenient to define the dimensionless correlation matrix corr^j(p) = coYij{p)/a{pi)a{pj). Geometrically, the 
correlation matrix can be envisioned as an ellipsoid in the 4-dimensional parameter space fulfilling x^ “ Xmin ~ 1* 
While the diagonal components of the matrix are unity by definition, the off-diagonal components are found between 
-I and I, and contain information about the correlation between different fitting parameters, which quantifies how 
well the function and the associated parameters are chosen to reproduce the data. A large degree of correlation results 
in a large redundancy between the parameters, a large fitting error and a small degree of confidence of the fitting 
procedure. We find that the bi-exponential function models accurately the data judging from the uncertainties being 
a few percent or less, see Sec. , as well as from the low degree of correlation between the parameters that was checked 
for every fit. 
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FIG. 10: (a) Sketch of the interface-fluctuation GaAs quantum dot embedded in an AlGaAs matrix (not to scale). The 

two-dimensional electron-hole pair (exciton) is coherently spread over the spatial extent of the quantum dot (green area). 
Exciton enhancement is achieved within the plane (grey arrows), while out of plane cooperative effects are destroyed by the 
close proximity of the GaAs-AlGaAs potential barrier (black arrows), (b) Band diagram in the transverse direction and the 
corresponding wave functions of the investigated exciton. The material parameters used in the calculation are taken from 
Ref. 50. (c) Plot of the in-plane exciton density |Tx(iCe, 0, 0)|^. 


VI. EXTRACTING THE SIZE OF THE EXCITON WAVE FUNCTION 


The presented experimental findings provide insightful information not only about macroscopic properties such as 
the oscillator strength and quantum efficiency, but also about the nanoscopic structure of the quantum-mechanical 
wave functions of the quantum dot. The out-of-plane uncorrelated electron and hole wave functions are computed 
with a tunneling resonance technique [54] and are plotted in Fig. 10(b) for the investigated quantum dot. While the 
microscopic structure of the out-of-plane wave functions can be computed because the number of atomic layers in 
the quantum well is known precisely (see Sec. ), the in-plane geometry is generally unknown because the quantum- 
well thickness fluctuations are spatially random. This nanoscopic information is then inferred from the superradiant 
enhancement of spontaneous emission S', where it can be shown (see Sec. ) that the quantum dot radius L is related 
to S via 


j- _ ciQw 

V2 \{i^h\i^e)\’ 


(28) 


From the measured value S ^ 4.3 an in-plane diameter 2L 24 nm is obtained. The resulting wave function 
O 5 0) is plotted in Fig. 10(c), where a strong correlation between the electron and hole position within the 
quantum dot is observed, which gives rise to superradiant emission. Our results emphasize that optical spectroscopy 
is a robust, non-invasive way of acquiring profound insight into the nanoscopic wave functions of quantum emitters. 
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FIG. 11: Decomposition of the in-plane exciton dynamics into (a) a center-of-mass and (b) a relative motion, (a) The center- 
of-mass dynamics describes the motion of the exciton center of mass in a two-dimensional harmonic potential with equally 
spaced eigenstates, cf. Eq. (29). (b) The electron-hole electrostatic attraction is captured by the relative-motion dynamics, 
which is analogous to the two-dimensional Hydrogen with a hyperbolic confinement potential. 


VII. THE IMPACT OF THERMAL EFFECTS ON SINGLE-PHOTON SUPERRADIANCE 

The central requirement for SPS is quantum dots larger than the exciton radius. The quantization energy AEqd 
scales inversely proportional to the quantum-dot size squared, AF^qd cx: 1 /“^, and, thus, decreases dramatically for 
large quantum dots. If AF^qd is comparable to the thermal energy /cpT, excited states of the exciton manifold become 
populated with a potentially detrimental effect on the superradiant enhancement of light-matter interaction. This 
can happen even at cryogenic temperatures owing to the large size of monolayer-fluctuation quantum dots. In the 
following, we estimate the impact of thermal effects and show that the temperature limits the oscillator strength that 
can be harvested from a quantum dot. In particular, oscillator strengths larger than about 100 cannot be resolved at 
the base temperature of our cryostat (7K). 

The process of thermal activation of excited states is intimately connected to the quantum-dot energy structure. We 
model the monolayer-fluctuation quantum dots as cylindrically symmetric with in-plane parabolic and out-of-plane 
step-like quantum confinement as explained in Sec. . The out-of-plane confinement is strong and thermal effects play 
no role at cryogenic temperatures. In the following, we address the weak confinement within the plane. Since the 
exciton dynamics can be decomposed into a center-of-mass (CM) and a relative (r) motion, cf. Eqs. (6-8), the exciton 
manifold contains contributions from both. The confinement potential is cylindrically symmetric and the ground state 
can therefore be denoted as |1s)^ while an excited state can comprise any other combination of the subsets, 

e.g., |2 s)qj^ |2s)qj^ |2p)r: etc. The relative motion is mathematically equivalent to the two-dimensional Hydrogen 

problem [26, 55] and is governed by the mutual electrostatic attraction between the electron and the hole, see Fig. 11, 
which is why it does not depend on the quantum-dot size. In this subspace, the relevant energy difference AE^ 
between the ground |ls)^ and first excited |2s)^ states equals roughly twice the excitonic Rydberg energy (not four 
times as for an ideal two-dimensional system, cf. the discussion around Eqs. (6-8)) and amounts to about 8meV. At 
cryogenic temperatures, thermal energies are much smaller (below ImeV) and thermal population of excited states 
can be safely neglected. We therefore conclude that only the ground state of the relative-motion dynamics |ls)^ is 
populated, and that thermal effects only affect the spatially extended center-of-mass dynamics. 

The center-of-mass motion is described by a particle (exciton) in a two-dimensional harmonic potential Vcm (R) = 
(1/2)MU^R^, cf. Fig. 11(a), where M is the exciton mass and the spring constant U is related to the quantum-dot 
size L via [29] U = EhjML?. The resulting energy eigenstates are given by [55] 

F;,z = (2n-|/|-l)fiO, (29) 

where n = 1, 2, 3,... and / = 0,±l,...,±(n — 1). The dipole selection rules dictate that states with I = 0 are bright 
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FIG. 12: Plot of the normalized transition strength JV(T) versus normalized thermal energy k^T/Ml for different numbers 
of bound states in a monolayer-fluctuation quantum dot. J\f starts decreasing monotonically at around hTt ~ 4/cbF, which 
decreases light-matter interaction. 


(superradiant) and all others are dark (subradiant). The present paper studies spontaneous emission from the ground 
state = |lnOz): which is the only non-degenerate superradiant state of this manifold. For thermal energies 

k^T comparable to AF^cm = excited states become populated and the relevant figure of merit for light-matter 
interaction is no longer the oscillator strength /, but rather the transition strength F because spontaneous emission no 
longer involves a two-level system [19]. The transition strength is related to the oscillator strength via a temperature- 
dependent factor A/'(T), i.e., F(T) = Af{T) x /, wherein Af{T) can be regarded as a normalized transition strength 
and quantifies the distribution of the excitonic population within the center-of-mass subspace. It was shown in Ref. 19 
that for a single excited state, |e), A/" is given by 


Ar{T) = 


1 + ^B(T) 

1 + B{T) ’ 


(30) 


where B(T) = exp (—AF^cm/^b^) is the Boltzmann factor and /e is the oscillator strength of the excited state. If 
|e) is dark, /g ~ 0, and the temperature is sufficiently high, S ~ 1, the transition strength is half of /, i.e., M = 1/2. 
Here we are interested in finding the threshold temperature at which a certain oscillator strength / of a monolayer- 
fluctuation quantum dot can still be resolved. This corresponds to the regime in which no excited states are populated 
and F ^ f. Since monolayer-fluctuation QDs have, in general, a large center-of-mass subspace, we generalize Eq. (30) 
with the help of Eq. (29) and by noting that the oscillator strength of s-type states equals / while all other states are 
dark 


Af{T) 


radiative contributions 
all contributions 




^n-1 ^2(n-l)- 

2^n=l 2^Z=-(n-l) ^ 


(31) 


where N denotes the number of center-of-mass states. This expression can be evaluated analytically yielding 


Af{T) = coth 



tanh 



(32) 


and is plotted in Eig. 12. Eor confinement energies hCl smaller than about 4 /cbT, excited states play a negligible role 
and F = f. This is the regime in which the oscillator strength can be reliably measured. Our experiments are in 
this regime because we observe negligible emission from excited states below saturation. We employ this criterion 
hfl = 4/cbT to estimate the maximum oscillator strength max resolved at a temperature T and obtain 

(here and in the following we consider the out-of-plane overlap ~ 1 for simplicity) 

fSym _ 4FEp 1 . X 

/max,th Mujalk^T' ^ ^ 

which leads to an oscillator strength of 170 at the base temperature of our cryostat of 7K. This calculation has been 
performed for an in-plane cylindrically symmetric quantum dot. We evaluate the same expression for a more realistic 
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FIG. 13: Plot of the transition strength F versus quantum-dot radius L at different temperatures. For small quantum dots, F 
increases quadratically with L due to superradiance. The transition strength saturates when the thermal energy is larger than 
the quantum-confinement energy, and becomes independent of quantum-dot size. 


asymmetric shape with an aspect ratio of 1:^ with ^ > 1 and find that the maximum oscillator strength is decreased 

bye 


/i 


max,th 


nsym 
J max,th 


(34) 


which is equivalent to Eq. (3) in the main text. We therefore conclude that oscillator strengths larger than about 
100 are unlikely to be resolved at the experimental conditions of the present work, and lower temperatures are 
required for pushing this limit. As demonstrated in Ref. 29, oscillator strengths of 1500 are realistically achievable 
in monolayer-fluctuation quantum dots with a radius of about 60 nm, where the electric-dipole approximation is 
reasonably good. Temperatures below 0.8 K would, however, be required to resolve them. For completeness, we note 
that lim^^oo M{T) = 1/N^ i.e., the transition strength is inversely proportional to the number of bound center-of-mass 
states in the quantum dot at high temperatures. 

A natural question arises regarding the transition strength as a function of quantum-dot size L. There are two 
competing processes that scale in opposite directions. On the one hand, the coherent enhancement of light-matter 
interaction (i.e., SPS) increases with cf. Eq. (14). On the other hand, the normalized transition strength M 
decreases monotonically with ex as seen from Eig. 12. Erom these heuristic arguments it follows that the 

size dependence cancels out and F saturates for a sufficiently large L, which is confirmed in the following quantitative 
analysis. We assume that the quantum dot has a large center-of-mass subspace A, i.e., coth (fillA/2/cpT) 1, and, 
from F = A/, obtain 


F = U 


(Al' 

yCLQW ^ 


tanh 


(- 


2h^ 




\MkBTL2 J 


(35) 


which is plotted in Fig. 13. The transition strength saturates for large L as expected from the aforementioned 
arguments, and reaches a maximum value of 

L^oo MujaQ /cpT 

which is independent of L and, interestingly, happens to equal Note that for very large L > 100 nm, 

deviations from the electric-dipole approximation, which are not accounted for in this study, further reduce the 
transition strength [29]. 


(36) 


VIII. RESULTS OF ALL MEASURED QUANTUM DOTS 

The measurement results for all quantum dots that we studied are presented in Tab. 1. They represent all measured 
quantum dots (9 in total), i.e., no data have been discarded after analysis. Eigure 14 shows the quantum dot 
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FIG. 14: Data of QD3 exhibiting the highest measured oscillator strength in this study, (a) Normalized spectra under heavy- 
hole (green) and 2s-shell excitation (blue) taken at half of the saturation power, (b) Associated decay curves from which we 
extract an oscillator strength of 96. A good fit to the data obtained with C-excitation (2s-excitation) is obtained with a double 
(single) exponential function in agreement with the theoretical model. 


(QD3) exhibiting the fastest decay rate of 11.1 ns“^, which corresponds to an oscillator strength of 96. All studied 
quantum dots have a large oscillator strength with an average value of 76, which constitutes an average superradiant 
enhancement of 4.4 over the strong confinement limit. The average quantum efficiency is 95 %. The oscillator strength 
and quantum efficiency have been estimated with an accuracy of about 2 % as discussed in Sec. . 


Quantum dot 

Frad (nS“^) 

Fnrad (uS ) 

Fsf (ps“^) 

/ 

ri{%) 

QDl 

8.3±0.1 

0.130 ±0.003 

360 ±30 

72.0 ±0.8 

99 ±2 

QD2 

8.35 ±0.07 

0.406 ±0.008 

33 ±2 

72.0 ±0.6 

96 ±2 

QD3 

11.1T0.2 

1.3±0.2 

160 ±30 

96 ±2 

90 ±2 

QD4 

10.5T0.2 

0.414 ±0.003 

47±2 

90 ±2 

96 ±2 

QD5 

7.6 ±0.2 

0.93 ±0.04 

106 ±8 

66 ±2 

90 ±3 

QD6 

9.7±0.2 

0.34 ±0.02 

8±2 

83 ±1 

97±2 

QD7 

7.13 ±0.08 

0.300 ±0.007 

13 ±2 

61.5±0.7 

96 ±2 

QD8 

8.13 ±0.07 

0.374 ±0.007 

15.7±0.9 

70.1 ±0.6 

96 ±2 

QD9 

8.3±0.1 

0.454 ±0.007 

106 ±5 

72.0 ±0.9 

95 ±2 


TABLE I: Data extracted from time-resolved measurements on all measured quantum dots: radiative decay rate Frad, non- 
radiative decay rate Fnrad, spin-flip rate Fgf, oscillator strength /, and quantum efficiency rf. QDl: data presented in most of 
the main article and here, QD2: data of the PLE, QD3: largest oscillator strength. The errors have been evaluated for every 
quantity as discussed in Sec. . 










